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Gravitational and dilaton radiation from a relativistic membrane 
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Recent scenarios of the TeV-scale brane cosmology suggest a possibility of existence in the early 
universe of two-dimensional topological defects: relativistic membranes. Like cosmic strings, 
oscillating membranes could emit gravitational radiation contributing to a stochastic back- 
ground of gravitational waves. We calculate dilaton and gravitational radiation from a closed 
toroidal membrane excited along one homology cycle. The spectral-angular distributions of 
dilaton and gravitational radiation are obtained in a closed form in terms of Bessel's functions. 
The angular distributions are affected by oscillating factors due to an interference of radiation 
from different segments of the membrane. The dilaton radiation power is dominated by a 
few lower harmonics of the basic frequency, while the spectrum of the gravitational radiation 
contains also a substantial contribution from higher harmonics. The radiative lifetime of the 
membrane is determined by its tension and depends weakly on the ratio of two radii of the 
torus. Qualitatively it is equal to the ratio of the membrane area at the maximal extension to 
the gravitational radius of the membrane as a whole. 



I. INTRODUCTION 

Now it is well understood that topological defects 
which could have been created at the early stages of cos- 
mological expansion may leave their footprints in the 
spectrum of the background gravitational waves ^ ^ 
which can be detected in the future experiments. Un- 
til recently main attention was paid to the effects due 
to cosmic strings. Cosmic strings, presumably serving 
as seeds for the formation of galaxies, were thought to 
be produced at the phase transition of the grand uni- 
fication scale H |[. Gravitational, dilaton and axion 
radiation emitted by oscillating strings was thorough! 
studied in a number of papers (see e.g. ||, |6[ 0, |F 
and references therein). Although a recent analysis of 
the available astrophysical data seems to rule out the 
cosmic string scenario of the formation of structures in 
the Universe, these calculations may still be of interest 
in the context of the fundamental string/M theory. 

Recently it was argued that the mass scale of the su- 
perstring/M theory may be as low as 10 TeV 0, [ll[] 
at the expense of existence of some non-small extra di- 
mensions. If true, this may drastically modify the full 
scenario of cosmological evolution. Based on this as- 
sumption various cosmological models were envisaged in 
which the observed Universe is treated as a three-brane 
embedded in a higher-dimensional world. According to 
this brane cosmology the standard model particles are 
confined to the brane while gravity lives in the full space. 
Such scenarios appeal to reconsider the problem of topo- 
logical defects which should be associated rather with 



solitons of the fundamental string theory than with those 
of the field-theoretical models of grand unification. 

Among the string theory solitons which may live on a 
three-brane there are one-dimensional strings and two- 
dimensional membranes. In the zero-thickness limit the 
dynamics of solitonic strings is governed by the Nambu- 
Goto action, while the membranes are described by a 
similar three-dimensional geometrical action (see e.g. 
||l2 , 13 1). Relativistic membrane in the four-dimensional 
spacetime geometrically are similar to domain walls en- 
countered in the usual field-theoretical models. How- 
ever their formation and dynamics in the string/brane 
cosmology are likely to be different from dynamics of 
the domain walls in the standard cosmological scenario. 
Here we do not discuss details of creation and evolution 
of membranes in the brane cosmology but concentrate 
on the problem of radiation. 

Radiation produced by one-dimensional defects - the 
relativistic strings - is likely to be similar to radiation 
from the cosmic strings, apart from a different attribu- 
tion of scalar and two-form fields involved. Radiation 
from an oscillating two-dimensional relativistic mem- 
brane seems not to have been studied in detail before, al- 
though some earlier estimates of gravitational radiation 
emitted by the cosmic domain walls are worth to be men- 
tioned |l(| [l~7|] . Massless boson modes to which the mem- 
branes are coupled in the four-dimensional space-time 
include a dilaton, a graviton and a Ramond-Ramond 
three-form. The latter is non-dynamical in four dimen- 
sions, so oscillating membranes may emit dilatons and 
gravitons. 
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In the classical limit one can describe radiation from 
a vibrating membrane in a standard way using the re- 
tarded solutions of the corresponding wave equations. 
Unfortunately, contrary to the case of strings, a general 
solution of the (highly non-linear) membrane equation 
of motion can hardly be obtained in a closed form. To 
simplify the problem one can restrict oneself by some 
particular excitation mode. Here we study the dilaton 
and gravitational radiation emitted by a toroidal oscil- 
lating membrane excited along a homology cycle with 
the radius much smaller than the second radius of the 
torus. Under such an assumption the membrane equa- 
tion of motion can be easily solved, and a closed ana- 
lytic expression for the spectral power of radiation can 
be obtained. We analyse in detail the spectral-angular 
distributions with an emphasis on the effects due to two- 
dimensional nature of a membrane. Hopefully this sim- 
ple model correctly describes basic features of the mem- 
brane radiation in more general case, although radiation 
due to more complicated excitation modes is certainly 
worth to be considered in the future. 

The plan of the paper is as follows. In Sec. II we 
derive the linearized equations for the scalar and gravi- 
tational radiation from a membrane source. In Sec. Ill 
the solution of the equations of motion for a free toroidal 
membrane is obtained under an assumption of the rela- 
tive smallness of the radius of an excited homology cy- 
cle. Sec. IV is devoted to an analysis of the dilaton 
radiation; we derive the spectral-angular distribution in 
terms of the Bessel's functions and investigate the role 
of higher harmonics of the basic oscillation frequency. A 
simplified expression valid in the high-frequency range 
of the spectrum is obtained in terms of the MacDon- 
ald's functions. In Sec. V the gravitational radiation is 
studied along similar lines. The results are summarized 
and discussed in Sec. VI. Some useful integrals involving 
Bessel's functions used in the calculation are given in the 
Appendix. 



II. BASIC EQUATIONS 

Consider a relativistic membrane x 11 = z^(a a ), fi = 
0, 1, 2, 3, propagating in the four-dimensional space-time 
with the metric and interacting with the dilaton field 
4> (our signature convention is (H )). The corre- 
sponding action can be written as a sum 



S — S m -\- Sf 
of the purely geometrical membrane action [fL2l O) 



and the field term 



y/jd 3 a, 



(1) 



(2) 



(3) 



Here fj, is the membrane tension parameter of a dimen- 
sion of the energy per unit area, a is the dilaton coupling 
constant. In our normalization the dilaton field <p an d 
the coupling constant a both are dimensionless. The 
membrane world-volume is parameterized by the inter- 
nal coordinates a a = (t, er, p), a = 0, 1, 2. The induced 
metric on the world- volume reads 



(4) 



where 7 = det ^ a b ■ 

Variation of the action (|^) over the space-time coordi- 
nates z^ gives the equation of motion 
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M e ~°Vn° 6 ^(*%»')} -0, (5) 



where r ) ah is the inverse metric on the world-volume. 
Up to the dilaton factor, this is a three-dimensional co- 
variant D'Alembert equation. It is worth comparing 
it with the two-dimensional string equation of motion. 
In the latter case one has three gauge degrees of free- 
dom corresponding to two diffeomorphisms and a Weyl 
transformation on the world-sheet. Consequently, three 
independent components of the world-sheet metric can 
be chosen flat, so that the equation of motion becomes 
linear. A general solution can be written as a sum of 
two arbitrary functions of null coordinates on the world- 
sheet: right and left movers. This is no more possible in 
the case of a membrane for which the number of gauge 
degrees of freedom is again three (diffeomorphisms of 
the world- volume) , while the number of the independent 
world-volume metric components is six. Therefore the 
membrane equation of motion remains essentially non- 
linear in any gauge. To make it explicitly solvable one 
needs additional assumptions to be made about excita- 
tion modes. 

Variation of the full action (Q) with respect to the dila- 
ton (j> gives the four-dimensional D'Alembert equation 
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(6) 



with the source term 



j(x) = G^ae-"* [ 6 * {x Z± {a)) ^rd 3 a. (7) 



Note that the effective coupling of the dilaton to the 
membrane is proportional to the dilaton exponential. In 
absence of the dilaton background this factor is of course 
negligible in the linearized limit relevant to the radiation 
problem. 

Finally, varying the action (jl]) with respect to the 
space-time metric g^u, one obtains the Einstein equa- 
tions 



1 



(8) 
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with the source term consisting of the membrane stress- 
tensor 



where the source term is 



J 4 (*-*(*)) 



yfytfa, (9) 



and the dilaton term. In what follows we neglect the in- 
teraction of the dilaton with gravity and omit the second 
term. The effective coupling of the membrane to gravity 
contains again the dilaton exponential factor which may 
be neglected in absence of the dilaton background. 

The action (Q) is invariant under the reparameteriza- 
tion of the world- volume a a — > <j a (<7 6 ), so three gauge 
conditions may be imposed on the world-volume metric 
7 a h. We will treat the gravitational field perturbatively 
on the flat background 



(10) 



where rj^ is the Minkowski metric, and choose the pa- 
rameterization of the world-volume ensuring diagonality 
of the metric j a b' 



0. 



*% = o, 



0. 



(11) 



Here and below all the contractions are performed using 
the flat metric rj^, and the following abbreviations are 
used i M = d T z^, z' M = d^z^, z M = d p z^. The induced 
metric then takes the form 

z 2 



lab= | z' 2 

z 2 



= \W¥, (12) 



where z 2 = z^z^ etc. Note that i M is a timclike vector, 
while and z M are spacelike, so the signature of ^ a b is 

H (hence det7 a b is positive). In the diagonal gauge 

the flat-space membrane equation of motion (m can be 
presented as the continuity equation (see, e.g., [12]): 



dPg 
P =0, 



(13) 



Pit 



-.z^pz' 2 z 2 , 



dr da dp 

where P£ are generalized momenta obtained by varying 
the membrane action over dx^/ da a : 

1 

i 

7r 

i 

These equations are still highly non-linear. They can be 
solved, however, if one assumes certain symmetric con- 
figurations of the membrane. 



• 2 -2 

Z Z , 



(14) 



-z»z 2 z n . 



In absence of the gravitational background field (10), 
the dilaton equation (|^) can be considered in the 
Minkowski space-time. Denoting the flat wave operator 
by a box we have 



j(x) = fiaG / Vz 2 z' 2 z 2 8 4 (x - z(a)) d 3 a. (16) 



Finally, the linearized Einstein equations for the met- 
ric perturbations can be written in the Fock-de Dondcr 
gauge 



U^ v = 16ttGT^, 



(17) 



where as usual 



1>v> = h^, - -h^, 9^ = 0. (18) 

In the linearized theory the membrane stress-energy ten- 
sor simplifies to the following expression 



z»z v z' 2 z 2 + z'»z' v z 2 z 2 + z»z v z 2 z' 2 \ 



x<5 4 [x - z{a)) (z 2 z' 2 z 2 y 1/2 d 3 a{l9) 



Note that the integral in the right hand side contains one 
integration more as compared with the case of a string 
and two integrations more compared with the point par- 
ticle. The extended nature of the membrane source gives 
rise to the interference effects in the angular distribution 
of the radiation power. 



III. TOROIDAL OSCILLATING MEMBRANE 



Our choice of the membrane model is motivated as 
follows. In view of complexity of the general situation, 
we intend to calculate radiation coming from the low- 
est excitation mode of the closed membrane. First the 
topology of the membrane 2-dimcnsional surface has to 
be specified. The simplest spherical membrane does not 
radiate in the main (spherical) excitation mode. So we 
choose a toroidal membrane in which case excitations 
naturally produce a non-vanishing quadrupole moment. 
We also assume that only one homology cycle of the torus 
is excited, this considerably simplifies the equations of 
motion. 

Specifying z-direction of the coordinate system to be 
along the symmetry axis of the torus we arrive at the fol- 
lowing parameterization of the membrane world-volume 

z a = t, z 1 = [R + r(r) cos p] cos a, 

z 2 = [R + r(r) cos p] sin a, z 3 — r(r) sin p, (20) 

where R =const is the large radius and r(r) is the vari- 
able small radius. Two parameters on the world volume 



(15) 



ere [0,2?!-], pe[0,27r] 



(21) 
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measure the angular distance along large and small ho- 
mology cycles respectively, and r is a time coordinate on 
the world-volume equal to the global time in the chosen 
gauge. This ansatz satisfies the equations of motion ( |l3| ) 
only if at any moment r one has 



R 2 > r 2 . 



(22) 



In this approximation the world-volume element is equal 
to 



^=R^(l-r*), 



(23) 



substituting this into the equations of motion (|13|), one 
finds that under the assumption ( p2| ) the full system re- 
duces to a single equation 

rr + (1 - f 2 ) = 0. (24) 

This equation has the following solution 

r = Q- 1 sm(flT), (25) 

where fl is a constant, and it is assumed for simplicity 
that r(0) = 0. The parameter is the oscillation fre- 
quency which is equal to the inverse radius of the smaller 
circle ro of the torus at the moments of its maximal ex- 
tension. 

In our approximation each loop a =const belonging 
to the membrane moves independently. Actually the so- 
lution (plf ) satisfies exactly the equations of motion for 
an infinite cylindrical membrane. For a toroidal mem- 
brane this solution is only the main approximation valid 
if R > r . 

Note that the function r in ( |20| ) may take both posi- 
tive and negative values. Indeed, two points on the loop 
a =const which correspond to opposite ends of the di- 
ameter p = 0, 7r move to each other while r(r) is de- 
creased, merge when r = 0, and then exchange their 
positions while r(r) is negative. Therefore a half-period 
of r-oscillations corresponds to the full period of the true 
oscillations of the membrane. In other words, the period 
of the membrane motion is 7r/f2 and not 2tt/£1. 

In our units the inverse angular frequency is equal to 
the small radius of the torus ro = fi -1 at the moments 
of its maximal extension. Therefore the quantity 



v = nR = —, 



(26) 



equal to the ratio of two radii of the toroidal membrane 
at the maximal extension, may be regarded as the shape 
parameter of the membrane. 



(|16|). Since the membrane motion is periodic with the 
period 7r/S7, the radiation power can be presented as a 
sum over the even harmonics of the frequency 



(27) 



The angular distribution of radiation with the frequency 
uj n is given by 



dP n ^_ , 

dn 2ttG 



^ |j(k,^)| 2 



(28) 



where k is the wave vector, k 2 ~ cj 2 , and the Fourier- 
-transform of the source current is defined as follows 



7T 

j(k, uj n ) = - Jdt Je^-^jiv, t)d 3 x. 



(29) 



Parameterizing the wave vector k by the spherical angles 
0,tp: 



k = uj n (sin 6 cos ip, sin 9 sin ip, cos ( 



(30) 



and, using in (|29|) explicit expressions ( |16| ) and (|20|) , one 
obtains after integration over the space-time variables 
t r: 



^ J d X I dp I dae 2n ^ sin 2 x , (31) 



where 

tp = [x ~ sm P cos ^ snl X ~ v cos(er — (p) sin 0) . (32) 

Note that the azimuthal angle (p enters into this expres- 
sion (as well as to all other formulas below) only in the 
combination a — ip, and thus <p can be eliminated by 
the shift of the integration variable a — > a — p. This 
means that all quantities are actually ^-independent, as 
it could be expected in view of the axial symmetry of the 
membrane. Assuming without loss of generality <p = 0, 
one has 



ip = (x ~ sm P cos $ snl X ~ v cos a sm ' 



(33) 



The integration in (|3^) can be performed in terms of 
the Bessel's functions as follows. First one integrates 
over x using the integral representation (|7C|): 



IV. DILATON RADIATION 

The computation of the scalar radiation starting with 
the D'Alembert equation ( |l5| ) is standard and amounts 
to the evaluation of the Fourier-transform of the current 



^2ni(x— z sin x) 



d X — J cos [2n (\ — z sin X )] d\ 



+i j sin [2n (x — zsmx)] dx = nJ2 n (2nz) + il, (34) 
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where z = sin pcosO. It can shown that the contribution 
proportional to the imaginary part I vanishes after inte- 
gration over p. Indeed, a constant shift of the integration 
variable \ X + 7r /2 leads to the following integral 



r/2 



d\ / dp sin [2m (x — z cos x)] 



(35) 



-ir/2 



Here the integrand changes a sign under reflection of the 
domain of integration \ — * ~Xi P ~~ * — P- Writing in ( p2] ) 
sin x = (1 — cos2x)/2 one can show that the il term 
gives no contribution in presence of cos2x as well. 

The second integration over p is performed using the 
formula (^) of the Appendix. Finally, integrating over 
a and using the recurrent relations J72] ) one obtains 

i(k, cj„) = -2ttVg5 JoG,) (j2(a;)tan 2 + J* (a:)) , 

(36) 

where y — 2nv sin 9 and x — n cos 

In this expression the Bessel's functions J n (x) are typ- 
ical for radiation from relativistic sources like the syn- 
chrotron radiation from a relativistic charge. Another 
Bessel's function factor Jo(y) 2 has a different origin: it 
describes interference effects due to coherence of radi- 
ation emitted by different segments of the membrane. 
This factor reflects an extended nature of the membrane. 

One can see that the angular distribution of the to- 
tal power diverges in both directions along the symme- 
try axis 9 = 0, 7r. For sin (9 = 0, the interference factor 
Jo(0) = 1, so using an asymptotic formula (|7S| ) we obtain 
for n ^> 1: 



dPn(O) 

dn 



32 
¥ 



1/3 



Gp 2 a 2 R 2 



.-2/3 



The sum over n diverges, but one can show that the 
solid angle into which the high n radiation is peaked de- 
creases with growing n, so that the total radiation power 
remains finite. 

The angular distribution of the n-th harmonic is given 

by 



dP n 

dn 



8GTr 3 (panR) 2 J 2 (y) ( 



tan 2 9 J 2 + J'\ 



(38) 



where the argument of the Bessel's functions J n and J' n 
is x = n cos 9. This is shown at Figs. 1,2 for n = 1, 2 and 
v = 10. Angular distribution of the first harmonic has 
a maximum near the equatorial plane. Oscillations are 
due to the interference factor Jo(y). Recall that our as- 
sumption (^2|) is equivalent to v 3> 1, so for non-small 9 
one can replace Jo{y) by the asymptotic expression (77): 



Mv) 



2iTvn sin 9 



cos ( 2irvn sin 9 



(39) 



The distribution of the second harmonic exhibits max- 
ima displaced from the equatorial plane to some interme- 
diate region. It is modulated by the interference factor 
with doubled frequency. 

Performing the angular integration one obtains the 
distribution over the harmonic number as shown at the 
Fig. 3. One can see that about 88% of total energy loss 
coresponds to the first, and about 7% to the second har- 
monic, while the contribution from higher harmonics is 
relatively small. 

The high frequency tail may be described using ap- 
proximations for the Bessel's functions which follow from 
the asymptotic formulas (Fl],F4|) valid for large n 



J n (na) 
J' n {na) 



A'i 



/3 



sin# 
sin 2 9 



fn . 



sin J (9 



(40) 



where K 1 / 3 and K 2 /3 are MacDonald's functions. Since 
K u decrease exponentially for large values of the argu- 
ment ([76]) , it is clear that the high-frequency part of the 
radiation is concentrated in a narrow angular regions 
around 9 = 0, n 



sin 9 < n 1/<3 . 



(41) 



This beaming is a typical relativistic effect. The main 
contribution to the total power comes from the angular 
region in which the argument of the Bessel's function 
My) is lar S e 



y 



2flRn 1/3 > 1, 



(42) 



so the interference factor is rapidly oscillating. There- 
fore one can further simplify the angular distribution for 
higher harmonics using the asymptotic formula (^) and 
averaging over the interference ripples as follows 



4{y) 



1 

■ny' 



(43) 



As a result, the averaged over ripples angular distribu- 
tion of radiation for n 3> 1 will be given by the formula 



dP n AGp 2 a 2 Rn 



dn 

where 



(k 2 1/3 (z) + K 2 /3 (z))\ (44) 



-sm 



(45) 



Note, that this representation remains qualitatively good 
even for relatively small n, see Fig. 2. 

In view of a rapid convergence of the integrals from the 
MacDonald's functions, we can integrate this averaged 
distribution over angles. Passing to the argument of the 
MacDonald's functions as the new integration variable 
and using the relation 



sin 8 9d9 



-z dz, 



(46) 



G 



where an additional factor sin 9 comes from eK7, one ob- 
tains for n 3> 1 



Pn = 



32I 1 Gfi 2 a A R 1 



(47) 



where 



h = I ( K 2 1/3 + K 2 /3 ) dz = 4.33, (48) 



Therefore the contribution from high harmonics falls off 
as n~ 2 . 

One can also investigate the dependence of the radi- 
ated power over the shape parameter v (the ratio of two 
radii of the torus). Numerical probes for different values 
of v show a rather weak dependence of the total radiation 
power on this shape factor (see Fig. 4) 



V. GRAVITATIONAL RADIATION 

Gravitational radiation power can be computed along 
similar lines. In addition one has to distinguish between 
the polarisation states. The angular distribution of ra- 
diation at the frequency u n can be presented as a sum 
over two circular polarisations 



dP n Glo: 



dfl 



where 



T±(k, <>= £ ±r^(k, u n ) 



(49) 



(50) 



7T 

T""(k, w n ) = - / dtj e ^«*- kr )T^(r, t)d 3 x. (51) 



Here the standard chiral graviton projectors are used 



e ± =e ± e ± 

1 



(52) 



(12) 

where e/t ' are two linearly independent space-like unit 
four-vectors orthogonal to the wave four-vector — 
(w„,k) and its spatially reflected dual fc M = (u> n ,— k), 



e^k" = e^ 2) ^ = e«e^ = 0, 



-1. 



(53) 



Alternatively one could choose the linear polarisation 
states e® , e®, corresponding to real and imaginary parts 
of H): 



V2 



(54) 



or, explicitly, 



p© - _ ( P (i)p(i) _ P ( 2 ) P ( 2 )) 
e® - — fe (1) e (2) +e (2) e (1)> l 



V2 

In the Lorentz frame in which e 



(55) 



(1,2) 




two unit vec- 
tors orthogonal to k can be chosen along 9 and tp direc- 
tions: 

o(l) 



,(2) 



(cos#cos<y9, cos 6 sin ip, — sin (9), 
(— ship, cos ip, 0). 



(56) 



As in the dilaton case, the chiral projections of the 
Fourier-transform of the stress tensor depend only on 
the difference a — tp, so without loss of generality we can 
set tp = 0. After some simplifications one finds 



T±(k, u n 



2^VL 



7T — 7T 27T 



JdxJdpJ dae 2 "^$±, (57) 



where 



(cos (j cos 9 cos p — sin 9 sin p ± z sin a cos p)' ! 



(58) 

and ip is given again by the Eq.(|33|). 

After the integration over a, p both chiral amplitudes 
lead to the same result. Indeed, purely imaginary terms 
in $ ± are proportional either to sin 2<r or to sin 2p; taking 
into account an explicit form of the phase factor ip, one 
finds that both these terms do not contribute to the inte- 
gral because of the antisymmetry of the integrand in ap- 
propriate variables. In view of the relation ( p4| ) between 
chiral and linear polarisation projectors this means that 
the gravitational radiation in any direction has only one 
linear polarisation component e®,. By similar reasoning 
one can show that the product of the first two terms in 
(pq) also vanishes after an integration over p, so $ ± can 
be replaced by the sum of the squares of each of the three 
terms. After some rearrangements, we find the following 
equivalent expression suitable for further integration 



$1 = 



sin 



1 + cos 2 9 



-(1 — 3 cos2p) H cos2er(l + cos2p). 

(59) 

Note that, contrary to the dilaton case, this function 
is independent on Xi (*- e - on time). Therefore the in- 
tegration over x gives again the expression (|34|), and 
the imaginary part of the integral will vanish after the 
integration over p; this remains true for terms in ( |59| ) 
containing cos 2p. The subsequent integration over p is 
performed using the integral (p0|). The integration over 
a directly amounts to the integral representation for the 
Bessel's functions d7q), and finally one obtains 



T + (k,LJ n )-- 



(1 



■s 2 9)J 2 {y) (j'l 
+ sin 2 9J (y) ( J 2 (l + 3/ cos 2 9) - 3 J' 2 ) 



ir 2 pR 
2fl 



- tan 2 9 J' 



(60) 
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where the argument of the Bessel's functions of the or- 
der n, n± 1 is a; = ncosO. Now an interference of ra- 
diation produced by different parts of the membrane 
is accounted by two different Bessel's functions Jo(y) 
and Jziy); apparently this is related to the graviton 
spin. Since J2(0) = 0, the first term in (60) vanishes at 
6 = 0, 7r, so does the second term containing an explicit 
factor sin 2 8. Thus, contrary to the dilaton case, the am- 
plitude of the gravitational radiation strictly along the 
symmetry axis is zero. Still the main part of radiation 
is concentrated near this direction, see Figs. 5,6. 

The maximal contribution to the total radiation power 
comes from the first harmonic, but its relative contribu- 
tion is smaller than in the dilaton case (Fig. 7). The high 
frequency tail can be described in terms of the MacDon- 
ald's functions. For large y corresponding to the main 
contribution to the totat power one has Jo ~ J% (see 
(y7\), this can be used before averaging over the ripples. 
An averaged distribution in the high frequency region is 
then obtained as follows 



dP n 4 t i 2 Rn 



sin 7 6 (32^/3^) - Kl /3 (z)) 2 . (61) 



dfl 9ir 2 n 
After the angular integration one gets 

^ Sn 2 RI 2 1 



where 



3Kf 



/3 



K 2 

•"■2/3 



dz = 1.24. 



(62) 



(63) 



This is similar to the previous result for the dilaton (44). 

Numerical integration for different v shows that the 
total radiation power depends on the ratio of the torus 
radii even weaker than in the dilaton case (Fig. 8). 



VI. CONCLUDING REMARKS 

Let us summarize our results. We have shown that 
the dilaton and gravitational radiation from a vibrating 
relativistic membrane exhibits features typical for ex- 
tended relativistic sources: the radiation power contains 
high harmonics of the basic frequency beamed along the 
symmetry axis of the torus and its angular distribution 
is substantially modulated due to interference. The rel- 
ative contribution of the high frequency tail depends on 
the spin and it is more pronounced in the gravitational 
case. For high n the averaged over ripples radiation 
power falls off as P„ ~ 1/n 2 both in the dilaton and 
gravitational cases. The high frequency part of radiation 
can be conveniently described in terms of the MacDon- 
ald's functions like in the case of ultrarelativistic point 
particles. 

The total power of the radiation from a toroidal mem- 
brane can be presented in the form 



where A ~ Rrg is the membrane area at the maximal ex- 
tension, and the shape factor f(v) is a function of the ra- 
tio of two radii of the torus (actually we observed rather 
weak variation of f(v)). On dimensional grounds it can 
be expected that this formula will be valid for other ex- 
citation modes up to some factor, which we set to unity 
f(v) ~ 1 for a qualitative estimate. Then in terms of the 
total mass of the membrane M = [iA we get 



P 



GM 2 



so the radiative lifetime may be estimated as 



M 
~P 



A 
GM' 



(65) 



(66) 



This quantity is proportional to the ratio of the area at 
the maximal extension to the gravitational radius of the 
membrane as a whole (an inverse speed of light factor has 
to be inserted in the ordinary units). In our model the 
membrane area at the moments of a maximal extension 
is 4ir 2 R/tt (recall that fl is equal to an inverse radius 
of the small circle of the torus, and we assumed that 
R I/O), therefore the quantity [jlR/VL is proportional 
to the total mass. Thus the membrane lifetime 



fiR/fl 1 
~P G~Il 



(67) 



depends only on the membrane tension. In contrary, the 
lifetime the radiating string loop is proportional to the 
loop size. 

In the string case the radiative power in terms of a 
total mass has the order 



ps 



GM 2 

IF 



(68) 



where L is the string length. Correspondingly, the ra- 
diative lifetime is 



rad 



M 

pst 



L 2 

GM' 



(69) 



P = G^ 2 Af(v), 



(64) 



Actually our toroidal membrane has a shape of a 'thick' 
string, following this similarity one can identify L ~ 
R, A ~ Rtq. Then it is easy to see that the membrane 
lifetime is R/ro ^> 1 times shorter than that of a string 
of an equal length. 

It can be expected that gravitational radiation from 
membranes will depend on the membrane topology. For 
a toroidal membrane the radiation coming from the main 
excitation mode (as we have assumed here) in non-zero, 
while in the case of the spherical membrane the main 
sprerically symmetric excitation mode does not con- 
tribute to the gravitational radiation at all. 

This work was supported in part by the RFBR grant 
00-02-16306. 
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VII. APPENDIX 



In the main text we have used the following integral: 



Evaluating the integral in the representation of the 
Bessel's functions 



7T 

J cos2/ipJ2 U (2asm p)dp = 7r cos 7r/i J^-^ (a) J v +n(a) 



7T 

Jn(y) — — J cos (na — y sin a) da (70) In particular, in view of the Eq.(|72|) one has 



(80) 



for y — nx and n > 1 in the stationary phase approxi- 
mation one finds 



Mux) « v^EM Kl/3 Ml h{1 _ ,)3/ 2 j , (71) 

where K is the MacDonald's function. For the Bessel's 
functions one has the following reccurcnt relations 

Jn±x{nx) = -J n T J' n , (72) 

X 

while for the MacDonald's functions with non-integer v 
zK' v + vK v = -zK x _ v , (73) 



(note that K_ v = K u ). Differentiating fl7l| ) and using 
( [73| ) one obtains in the leading approximation for large 
n the corresponding formula for the derivative of the 
Bessel's function 

J' n {nx) a ^^2/3 (j^n{l - x)^ . (74) 

For small arguments the MacDonald's function behaves 
as 



K u {z) 
while for large ones 



-JL_r-i(i-„)(H 

2 sin ttv \z 



v W 2z 6 ' 



For Bessel's functions at large y one has 

Mv) ~ \lzr. cos ( 



2 / pn 7r 
— cos I y 
iry 



(75) 



(76) 



(77) 



Passing in the Eqs.( |7l] , |74] ) to the limit x — > 1, in which 
case the asymptotic behavior ( |75| ) holds, one obtains for 
n > 1 

/o \ 1/3 

J n (n) « 3~ 2 / 3 r- 1 (2/3) f -J , (78) 



J' tt (n) ^S-Var-^l/S) 



2/3 



7T 

y J2n(2nxsinp) cos2pdp = tt (j'^ n (nx) -J^(nx] 



(81) 
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FIG. 3: Distribution of the dilaton radiation over the FIG. 4: Total radiation power as a function of the 

harmonic number n. shape parameter v = flR. 
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FIG. 5: Angular distribution of the first harmonic FIG. 6: Angular distribution of the 4-th harmonic 

the gravitonal radiation (Pq = G[i 2 R/Q) the gravitational radiation and its averaged approx- 

imation by MacDonald's functions (see Eq. (fill), 
P = Gn 2 R/n. 




FIG. 7: Distribution of the gravitational radiation FIG. 8: Total gravitational radiation power as a 

over the harmonic number n. function of the shape parameter v = QR. 



